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Abstract — A  unified  treatment  of  several  least  squares 
(LS)  algorithms  is  presented  for  bearings-only  tracking  of  a 
target  moving  at  constant  acceleration.  The  close  link  be¬ 
tween  the  maximum  likelihood  (ML)  estimator  and  other 
nonlinear  and  “linearized”  LS  algorithms  is  explored  un¬ 
der  the  assumption  of  Gaussian  bearing  noise.  In  this  con¬ 
text,  a  new  asymptotically  unbiased  closed-form  instrumen¬ 
tal  variables  (IV)  algorithm  is  derived.  Reduced-bias  total 
least  squares  (TLS)  and  constrained  TLS  (CTLS)  algorithms 
are  developed.  The  equivalence  of  the  ML  algorithm  to  the 
structured  TLS  (STLS)  algorithm  is  established.  Simulation 
examples  are  provided  to  demonstrate  the  improved  perfor¬ 
mance  of  the  IV  and  TLS  estimators  vis-a-vis  the  pseudolin- 
ear  estimator. 

I.  Introduction 

The  paper  presents  a  unified  treatment  of  several  classical 
and  new  LS  algorithms  for  bearings-only  tracking  of  a  tar¬ 
get  moving  at  constant  acceleration.  Tracking  of  constant- 
velocity  targets  and  stationary  target  localization  are  spe¬ 
cial  cases  of  constant-acceleration  target  tracking.  The 
key  to  target  tracking  is  the  estimation  of  target  trajectory 
from  noisy  bearing  and  own-ship  location  measurements. 
For  targets  moving  at  constant  acceleration,  the  target  tra¬ 
jectory  is  defined  by  the  target  motion  parameters,  viz., 
position,  velocity  and  acceleration.  Target  tracking  has 
many  civilian  and  military  applications,  such  as  air  traffic 
control,  surveillance,  etc. 

Firstly  we  derive  the  classical  ML  estimator  for  the  tar¬ 
get  motion  parameters.  The  ML  estimator  becomes  a  non¬ 
linear  LS  estimator  for  zero-mean  Gaussian  distributed 
bearing  noise.  While  the  ML  estimator  is  optimal,  it  does 
not  have  a  closed-form  solution  because  of  the  nonlinear 
relationship  between  the  measurements  and  the  unknown 
parameters.  Based  on  the  classical  ML  estimator,  a  lin¬ 
earized  weighted  LS  estimator  is  derived  under  the  as¬ 
sumption  of  small  bearing  noise,  which  yields  a  closed- 
form  solution.  This  algorithm  is  simply  an  extension  of 
the  target  localization  algorithm  developed  by  Stansfield 
in  [1]  to  tracking  of  moving  targets.  A  pseudolinear  es¬ 
timator  is  obtained  from  the  Stansfield  estimator  by  re¬ 
moving  the  weighting  matrix.  An  alternative  derivation 
for  the  pseudolinear  estimator  based  on  orthogonal  vec¬ 


tor  representation  yields  the  equivalent  orthogonal  vectors 
(OV)  estimator.  The  linear  LS  estimators  are  known  to 
be  biased  [2],  A  new  closed-form  asymptotically  unbi¬ 
ased  IV  estimator  is  proposed  in  the  paper.  The  instru¬ 
mental  variables  are  obtained  from  another  estimate  such 
as  LS  or  TLS.  Unlike  the  algorithms  in  [2,  3],  no  recur¬ 
sive  computations  are  required  for  the  new  IV  estimator. 
TLS  and  CTLS  estimators  are  developed  to  improve  the 
performance  of  the  OV  estimator  by  attempting  to  miti¬ 
gate  errors  in  both  the  system  matrix  and  the  data  vec¬ 
tor.  A  structured  TLS  estimator  is  also  formulated  and  the 
equivalence  between  the  STLS  and  ML  estimators  is  es¬ 
tablished.  The  performance  improvement  achieved  by  the 
TLS  and  IV  estimators  compared  with  the  OV  estimator  is 
demonstrated  with  simulation  examples.  The  simulations 
also  show  the  capability  of  the  IV  estimator  to  outperform 
the  ML  estimator. 

II.  Problem  Formulation 

The  two-dimensional  target  tracking  problem  using  bear¬ 
ing  measurements  only  is  depicted  in  Fig.  1.  The  objective 
of  bearings-only  target  tracking  is  to  identify  the  target  lo¬ 
cation  pk  from  noisy  bearing  and  observer  position  mea¬ 
surements  over  a  finite  time  interval  0  <  k  <  N  —  1. 

In  Fig.  1,  the  relation  between  the  bearing  angle,  ob¬ 
server  position  and  target  location  is  given  by  the  follow¬ 
ing  nonlinear  equation 

Ok  =  tan-1  k  =  0,...,N  -1  (1) 

Axk 

where  A yk  —  PyA  Aa.^  Px,k  pk  — 

[Px,k,Py,k]T  is  the  target  location  vector  and  rk  =  [rx,kt 
Ty,fe]7  is  the  observer  position  at  k. 

We  make  the  following  assumptions  about  the  target 
tracking  problem: 

•  The  target  is  moving  at  constant  acceleration.  Let 
pQ  and  Vo  denote  the  target  position  and  velocity 
vector  at  k  =  0,  and  a  be  the  constant  target  ac¬ 
celeration  vector.  Assuming  that  the  bearing  and 
observer  position  measurements  are  taken  at  reg¬ 
ular  time  instants  tk  =  kT  where  T  is  the  sam¬ 
pling  interval,  the  target  location  at  time  tk,  k  = 
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Fig.  1.  Two-dimensional  bearings-only  target  tracking  ge¬ 
ometry. 


III.  LS  Tracking  Algorithms 
A.  ML  Estimator 

Assuming  that  the  observer  position  errors  are  zero,  that 
is,  wk  =  0  in  (4),  and  using  the  Gaussianity  assumption 
for  bearing  measurement  noise,  the  likelihood  function  for 
the  bearing  measurements  can  be  written  as 

=  (27r)JV/2|iT|i/2 

where 

9  =  [§o,  O2,  ■  ■  ■ ,  On-i]T 

is  the  A"  x  1  vector  of  noisy  bearing  measurements, 

o(0  =  [e0ld2,...,eN-1]T 


is  the  TV  x  1  vector  of  bearing  angles  which  are  dependent 
on  the  target  locations  pk  and  the  target  motion  parame- 
(2a)  ters  £  through  (1)  and  (2),  K  =  diag {a2no,. . . ,  a2nN_1) 
is  the  TV  x  TV  diagonal  covariance  matrix  of  the  bearing 

(2b) 

measurement  errors,  and  \K\  is  the  determinant  of  K. 

The  maximum  likelihood  estimator  of  the  target  mo¬ 
tion  parameters  £ML  is  obtained  from  maximization  of  the 
likelihood  function  p(0\£,)  over  all  possible  £  [5].  To  sim¬ 
plify  the  maximization  problem,  the  log-likelihood  func¬ 
tion  is  used: 


and 

>0 

£  =  v0 

a 

is  the  6x1  target  motion  parameter  vector  to  be  es¬ 
timated.  Given  an  estimate  of  £,  the  target  locations 
can  be  obtained  from  (2). 

•  The  bearing  measurements  are  subject  to  indepen¬ 
dent  zero-mean  Gaussian  noise: 

h  =  0k  +  nk,  nk  ~  A/'(0,  <r£j  (3) 

where  the  0k,  k  =  0, . . . ,  TV  —  1,  are  the  bearing 
measurements  available  for  target  tracking,  and  nk 
is  a  Gaussian  random  variable  with  zero  mean  and 
variance  a^k . 

•  The  observer  position  measurements  are  subject  to 
independent  bivariate  Gaussian  noise: 

rk  =  rk  +  wk,  wk  ~  Af(0,  C)  (4) 


lnp(0|£)  =  -^ln((27r)JV|iT|) 

-  5(0 -*(£))■ 

Noting  that  the  first  term  of  the  log-likelihood  function  is 
independent  of  the  maximization  of  the  log-likelihood 
function  can  be  achieved  by  minimizing  the  second  term 
only  after  sign  inversion,  leading  to 


£ml  =  arg  min  *7ml(£)  (5) 

£GR6 

where  Jml(£)  is  the  ML  cost  function  given  by 

Jml(€)  =  \{o  -  e{i))TK-\e  -  0(£)).  (6) 


The  minimization  of  Jml(£)  over  £  is  in  fact  a  nonlin¬ 
ear  LS  problem.  The  solution  to  this  minimization  prob¬ 
lem  satisfies 


<9Jml(£) 


(7) 


where  the  covariance  matrix  C  is  diagonal,  i.e.,  the  A  closed-form  solution  to  (7)  does  not  exist.  This  re- 
errors  in  x  and  y  coordinates  of  the  observer  posi-  quires  the  use  of  numerical  gradient-based  search  tech- 
tion  measurements  are  independent.  niques.  The  gradient-based  search  techniques  will  yield  a 

unique  solution  for  £ML  if  the  ML  cost  function  is  convex 
•  The  observer  trajectory  is  such  that  the  target  is  ob-  with  a  unique  global  minimum.  The  convexity  of  the  ML 
servable.  cost  function  was  stated  in  [6]  without  proof. 
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B.  Stcinsfield  Estimator 


Assuming  that  the  bearing  errors  are  very  small,  i.e.,  Or  — 
0k  ss  0,  we  can  write  [7] 

dk~ekKi  sin(0k  -  0k). 


Substituting  this  approximation  into  (6)  yields  the  follow¬ 
ing  cost  function 


JV-l 

Js(p)  =  E 

fc=0  n« 

1 


sin2 (Ok  -  Ok) 


=  l(FZ-b)TW-\FZ-b) 
where  F  is  the  N  x  6  matrix  defined  by 


F  = 


ao  Mq 
aj  Mi 


b  is  the  N  x  1  vector 

b  = 


ctk  — 


sin  0k 
—  cos  0k 


Nx6 


a-oro 

a^r  i 


L“Ar-irJV-iJ 

and  W  is  the  N  x  N  diagonal  weighting  vector 


W  = 


'dial  o 


2^2 


dfa; 


0 


N-l<JnN-i  J 


(8a) 

(8b) 


(9) 


(10) 


(ID 


Here  dk  =  yj{ Axk )2  +  (A yk)2  denotes  the  target  range 
from  the  observer  location  at  time  instant  k.  The  weight¬ 
ing  matrix  is  therefore  dependent  on  £. 

To  obtain  a  linear  LS  solution  to  the  bearings-only  tar¬ 
get  tracking  problem,  the  dependence  of  W  on  the  target 
motion  parameter  vector  £  is  ignored  by  using  the  true 
range  values  between  the  observer  and  the  target,  assum¬ 
ing  that  this  information  is  available.  This  simplification 
allows  us  to  treat  (8b)  as  a  weighted  LS  problem  with  a 
closed-form  solution  given  by 

£s  =  arg  min  Js(£)  (12a) 

£eR6 

=  (. FTW-1F)-1FTW~1b .  (12b) 


This  estimator  is  an  extension  of  the  Stansfield  target  lo¬ 
cation  estimator  [1]  to  moving  targets.  We  assume  that 
the  matrix  F  is  full-rank  as  required  by  the  observability 
assumption. 


C.  OV  Estimator 

In  (3)  we  formulated  the  target  tracking  problem  by  refer¬ 
ring  to  a  statistical  model  for  the  bearing  measurement  er¬ 
rors.  Alternatively  we  can  establish  an  orthogonal  vector 


n 

Observer 

Fig.  2.  Relationship  between  noisy  bearing  vector  sk  and 
the  target  location  pk. 


sum  relationship  between  the  true  and  measured  bearing 
vectors  as  shown  in  Fig.  2  where  sk  is  the  noisy  (mea¬ 
sured)  bearing  vector  emanating  from  the  observer  posi¬ 
tion  rk  and  makes  an  angle  of  0k  with  the  x-axis,  sk  is  the 
bearing  vector  between  rk  and  the  target  pk ,  and  ek  is  the 
orthogonal  error  vector  (note  that  el sk  =  0). 

From  Fig.  2  we  have 

pk  =  rk  +  sk  (13a) 

=  rk  +  sk  +  ek  (13b) 


where  the  orthogonal  vector  ek  is  defined  by 

ek  =  ||sfc||2sin(0fc  -  0k) 

=  dk  sin  rife  ak. 


sin  0k 
—  cos  0k 


(14) 

(15) 


Here  dk  =  ||  sk  ||  2  is  the  range  between  the  observer  vector 
Tfe  and  the  target  pk  defined  after  (11),  nk  =  0k  —  0k  is 
the  Gaussian  bearing  noise  in  (3),  and  ak  is  the  unit  vector 
orthogonal  to  sk  defined  in  (9). 

To  eliminate  the  noisy  bearing  vector  sk  from  the  equa¬ 
tion,  multiply  (13b)  through  with  the  transpose  of  the  or¬ 
thogonal  unit  vector  ak,  yielding 

akPk  =  akrk  +  m  (16) 

where  r]k  =  dk  sin  nk  is  a  zero-mean  nonlinear  Gaussian 
noise.  Concatenating  (16)  for  k  =  0, . . .  ,N  —  1,  we  get 

F£  =  b  +  i 7  (17) 

where  F  and  b  were  defined  in  (9)  and  (10),  respectively, 
and  r)  =  [770, ... ,  rjN-i]T ■  An  LS  solution  to  (17)  is  given 
by 

€ls  =  argmin  || F£  -  6|| \  (18a) 

CeR6 

=  ( FTF)~1FTb  (18b) 

which  is  referred  to  as  the  orthogonal  vectors  (OV)  esti¬ 
mator  or  the  pseudolinear  estimator  [2] . 

Under  the  small  bearing  error  assumption  (i.e.,  nk  ~ 
0),  the  errors  r]k  can  be  approximated  by  rjk  «  dknk 
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which  suggests  a  Gaussian  rjk  with  variance  d\o/lk.  To 
handle  errors  rjk  with  different  variances,  a  weighted  LS 
solution  can  be  constructed: 

argmin  \\W~1^2(F£,  —  fo) || | 

£GR6 

which  is  identical  to  the  Stansfield  solution  in  (12). 

The  diagonal  weighting  matrix  W  aims  to  improve 
the  LS  solution  by  adjusting  the  contributions  of  equations 
concatenated  in  (17)  according  to  their  noise  variances. 
The  diagonal  entries  of  W^1^2  (i.e.,  the  weights)  are  de¬ 
termined  from  the  reciprocal  of  noise  standard  deviations. 
If  the  range  information  dk  is  not  available,  the  OV  esti¬ 
mator  is  preferred  over  the  Stansfield  estimator  as  it  does 
not  require  the  use  of  the  weighting  matrix.  In  practical 
applications,  the  range  information  for  the  target  is  not 
usually  available.  In  the  absence  of  this,  the  weighting 
matrix  is  typically  replaced  by  an  identity  matrix  (assum¬ 
ing  identical  bearing  noise  variances),  which  reduces  the 
Stansfield  estimator  to  the  OV  estimator. 

D.  IV  Estimator 

The  mean  of  the  OV  estimate  is 

E{ZLS}  =  Z-E{(FTF)-1FTr1} 

where  —E{(FTF)~1FTrj}  is  the  estimation  bias.  Even 
though  we  have  E{rj}  =  0,  the  estimation  bias  is  gener¬ 
ally  non- zero  because  F  and  77  are  correlated  (the  bearing 
noise  rik  appears  in  both  F  and  77)  [8].  The  OV  solution 
is  obtained  from  the  normal  equations 

FtF£ls  =  FTb. 

Since  the  reason  for  the  bias  is  the  correlation  between 
F  and  77,  we  can  consider  the  following  modified  normal 
equations: 

GtF£w  =  GTb 

where  G  is  the  matrix  of  instrumental  variables.  If  G  is 
chosen  such  that  E{GT  F}  is  nonsingular  and  E{GTrj }  = 
0,  then  £lv  will  be  asymptotically  unbiased,  i.e.,  /i’{£LS }  = 

£  as  N  — >  00. 

A  simple  and  appealing  choice  for  the  instrumental 
variable  matrix  G  is  first  to  obtain  an  estimate  for  the  tar¬ 
get  motion  parameters  using  one  of  the  LS  algorithms.  In 
view  of  its  improved  estimation  performance,  we  favour 
the  TLS  solution  £TLS  which  is  discussed  in  the  next  sec¬ 
tion.  The  target  location  estimates  pk  =  [px,k,Py,k]T  are 
then  obtained  from  (2),  using  £-ns.  Given  the  pk,  we  can 
construct  the  instrumental  variables  matrix: 


0-0  Mo/ do^no 

sin  9  k 

1 

1 

s 

<M  g 

7 

<r,b 

7 

1 

<s 

_ 1 

1  & k  — 

—  cos  9k 

where  dk  =  || pk  -  rk ||2  and  9k  =  tan  1  The 

asymptotically  unbiased  IV  estimator  is  then  given  by 

&v  =  ( GTF)~1GTb . 


E.  TLS  Estimator 

As  common  with  all  LS  solutions,  the  OV  estimator  im¬ 
plicitly  assumes  that  only  b  is  subject  to  error.  In  fact,  be¬ 
cause  of  the  presence  of  errors  in  bearing  measurements, 
the  system  matrix  F  is  also  subject  to  error.  To  improve 
the  accuracy  of  the  OV  estimator,  the  concept  of  total  least 
squares  (TLS)  can  be  invoked  to  mitigate  errors  in  both  F 
and  b  [9].  Formally,  the  TLS  estimate  £-ns  is  given  by  the 
solution  of  the  following  constrained  optimization  prob¬ 
lem  [10,  11]: 


min  ||£[A,-5]T||f  (19) 

(F+A)|TLS=b+5 

where  L  =  diag(7i, . . .  ,£jv)  is  an  N  x  N  diagonal  weight¬ 
ing  matrix,  T  =  diag(£i, . . .  ,£7)  is  a  7  x  7  diagonal 
weighting  matrix,  and  ||  •  \\p  denotes  the  Frobenius  norm. 
According  to  the  optimization  problem  in  (19),  the  TLS 
solution  is  obtained  by  adding  minimal  perturbations  to  F 
and  b  so  that  the  perturbed  matrix  equation  is  consistent. 

The  TLS  solution  can  be  obtained  from  a  singular  value 
decomposition  (SVD)  of  the  augmented  N  x  7  matrix 
L[F,  - b\T : 

L[F,-b]T=  UT,Vt  (20a) 

7 

=  ^  (20b) 


where  U  =  [ukl . . . ,  7x7]  is  an  N  x  7  unitary  matrix,  S  = 
diag((7i, . . . ,  >77)  is  the  7x7  diagonal  matrix  of  ordered 
singular  values  (i.e.,  o\  >  •  •  •  >  (77  >0),  and  V  = 
[77 1 , ,  777]  is  a  7  x  7  unitary  matrix.  The  TLS  estimate 
is  given  by 


£tls  — 


£  7^77 


tlUl7 


£6^67 


(21) 


where  v~j  =  [77-17,  •  •  • ,  is  the  seventh  column  of  V. 
Equation  (21)  assumes  that  the  smallest  singular  value  is 
unique,  i.e.,  >  <77.  This  assumption  may  be  violated 

on  rare  occasions.  If  this  happens,  alternative  solutions 
can  be  obtained  (see  [11]  for  details). 

The  TLS  solution  in  general  exhibits  smaller  estima¬ 
tion  bias  than  the  LS  solution  [12].  In  the  absence  of  bear¬ 
ing  and  observer  location  errors,  we  get  F0£  =  bQ  where 
F o  and  bQ  are  obtained  from  true  bearing  angles  and  ob¬ 
server  locations.  The  noisy  matrix  F  and  vector  b  are  re¬ 
lated  to  Fa  and  bQ  through  F  =  Fa  +  A0  and  b  =  b0  +  d0 
where,  for  77.7.  «  0, 


UqMq 

cos  9k 
sin  9k 

A  0  « 

5 

1 

v oro 


Fn-^vn-  1 
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The  weighting  matrices  can  be  chosen  by  taking  into  ac¬ 
count  the  error  variances  in  F  and  b.  If  we  set  L  =  I  and 
T  =  diag(l,  1,  e),  where  e  >  0,  and  let  e  — >  0,  the  TLS 
solution  converges  to  the  LS  solution. 

F.  CTLS  Estimator 

The  perturbed  system  matrix  F  +  A  and  data  vector  b  + 
8  that  result  from  TLS  do  not  retain  the  structure  of  the 
original  system  matrix  and  data  vector  defined  in  (9)  and 
(10).  The  accuracy  of  the  TLS  estimate  can  be  improved 
if  the  TLS  perturbations  are  forced  to  obey  the  structure  of 
the  system  matrix  and  data  vector.  The  imposition  of  these 
constraints  on  the  TLS  solution  leads  to  a  constrained  TLS 
(CTLS)  problem  that  can  be  formulated  as 


x-axis  (km) 


min||i[A,  <5]T’||  (22a)  Fig.  3.  Simulated  target  tracking  geometry. 


subject  to 

(F  +  A)£CTLS  =  b  +  8  (22b) 

o,q  0  0 

df  tidj  4-af 

(22c) 

+2 

||&i||2  =  l,  i  =  0,  ■  ■  ■  ,N  —  1  (22d) 

d'o't'o 

b  + 8=  :  (22e) 


subject  to 

F(0  +  8)£stls  =  b{9  +  5)  (23b) 

where  L  is  an  N  x  N  diagonal  weighting  matrix  and 

aT(cj)0)M0 

F(4>)  =  \  (24a) 

_aT  {4>n- l )  Mn- i_ 
aT(</>0)r0 

b(4>)  =  :  (24b) 

_aT(0Ar_i)rAr_i_ 


F  +  A  = 


A  detailed  treatment  of  CTLS  can  be  found  in  [9]. 

A  closed-form  solution  to  (22)  is  not  available  because 
of  the  nonlinearity  of  constraints.  We  therefore  seek  an 
iterative  numerical  solution  based  on  nonlinear  program¬ 
ming  and  the  method  of  successive  projections  [9].  The 
idea  behind  successive  projections  is  to  obtain  the  TLS  so¬ 
lution  and  then  to  project  it  to  vectors  satisfying  the  given 
constraints  using  a  least-squares  criterion,  and  repeat  the 
whole  procedure  until  TLS  converges  to  a  solution  satis¬ 
fying  all  the  constraints. 

The  method  of  successive  projections  was  shown  to 
produce  a  vector  sequence  that  always  contains  a  subse¬ 
quence  that  converges  to  a  vector  satisfying  the  constraints 
of  the  optimization  problem  under  some  mild  conditions 
on  the  projections  [13]. 

G.  STLS  Estimator 

The  CTLS  algorithm  attempts  to  maintain  the  structure 
of  the  system  matrix  and  the  data  vector  while  perturbing 
them  to  achieve  consistency.  The  constrained  optimiza¬ 
tion  problem  can  be  recast  as 

min  1 1 L8 1|  \  (23a) 

s 


with  <p=  [(j>0, . . .  ,4>n-i]T  anda(^)  =  [sin0i;  -  cos^]T 

Equation  (24)  makes  explicit  the  nonlinear  parameteriza¬ 
tion  of  F  and  b  in  terms  of  N  bearing  angles  0.  The 
constrained  optimization  problem  in  (23)  is  known  as  the 
structured  total  least  squares  problem,  and  its  solution  has 
been  addressed  in  [  14,  15]. 

The  STLS  problem  is  in  fact  identical  to  the  ML  prob¬ 
lem.  The  equivalence  between  the  STLS  and  ML  solu¬ 
tions  can  be  seen  by  setting  L  =  K  1  (2,  which  reduces 
(23a)  to  the  minimization  of  the  ML  cost  function  in  (6), 
and  by  rewriting  (1)  as 

sin  9k  _  py,k  -  ry,k 
cos  9k  px,k  —  rx>k 

which  shows  that  the  ML  solution  satisfies  the  consistency 
constraint  in  (23b). 

IV.  Simulation  Examples 

In  the  simulation  examples,  we  use  the  target  tracking  ge¬ 
ometry  shown  in  Fig.  3.  The  observer  trajectory  consists 
of  two  constant  velocity  legs.  The  target  moves  at  a  con¬ 
stant  acceleration  with  motion  parameters  pQ  =  [50, 100]T, 
Vq  —  [—  2,  —  8]t  and  a  =  [— 0.5, 1.5]T.  The  bearing 
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y-axis  (km) 


x-axis  (km) 

Fig.  4.  Average  of  target  track  estimates. 


TABLE  I  MSE  ( <jw  =  0.01) 


OV 

TLS 

IV 

ML 

0.1° 

4.8 

3.5 

3.0 

3.0 

0.3° 

122.0 

33.4 

28.1 

28.2 

0.5° 

554.9 

99.1 

79.2 

79.8 

O 

o 

1267.1 

191.5 

157.6 

159.6 

0.9° 

2128.5 

332.8 

294.4 

302.7 

TABLE  II  Bias  of  p0  Estimate  (aw  =  0.01) 

OV 

TLS 

IV 

ML 

0.1° 

1.18 

0.17 

0.01 

0.01 

0.3° 

9.64 

1.44 

0.11 

0.16 

0.5° 

22.08 

3.67 

0.38 

0.55 

O 

o 

34.17 

5.72 

1.21 

1.46 

0.9° 

44.77 

8.97 

1.21 

1.60 

measurements  are  taken  at  tk  =  kT,  where  T  =  0.5 
and  k  =  0, . . . ,  N  —  1.  For  the  moving  target,  N  = 
40  bearing  measurements  are  collected  by  the  observer  at 
marked  locations  in  Fig.  3.  The  bearing  measurements 
are  subject  to  i.i.d.  zero-mean  Gaussian  noise  with  vari¬ 
ance  cr^.  The  observer  location  measurements  are  cor¬ 
rupted  by  i.i.d.  zero-mean  bivariate  Gaussian  noise  with 
covariance  matrix  C  =  diag(<r^,,  <rj).  The  target  motion 
parameters  were  estimated  using  the  OV,  TLS,  IV  and  ML 
estimators.  The  results  of  1000  Monte  Carlo  simulations 
for  mean-squared  error  (MSE)  and  initial  velocity  estima¬ 
tion  bias  are  listed  in  Tables  I  and  II  for  different  bearing 
noise  standard  deviations  an  with  the  observer  position 
error  standard  deviation  fixed  at  <jw  =  0.01km.  Fig.  4 
shows  the  average  of  target  track  estimates  for  different  al¬ 
gorithms  at  <jn  =  0.5°.  For  the  target-observer  encounter 
depicted  in  Fig.  3,  the  IV  estimator  yields  the  best  accu¬ 
racy  narrowly  outperforming  the  ML  estimator.  The  TLS 
estimator  appears  to  have  a  much  improved  performance 
compared  with  the  OV  estimator. 
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